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. . Current rocket technology can not send the spaceship very far, because the amount of the chemical 

G\ \ fuel it can take is limited. We try to use dark matter (DM) as fuel to solve this problem. In this 

work, we give an example of DM engine using dark matter annihilation products as propulsion. The 
acceleration is proportional to the velocity, which makes the velocity increase exponentially with 
time in non-relativistic region. The important points for the acceleration are how dense is the DM 
density and how large is the saturation region. The parameters of the spaceship may also have 
great influence on the results. We show that the (sub)halos can accelerate the spaceship to velocity 



10 5 c ~ 10 3 c. Moreover, in case there is a central black hole in the halo, like the galactic center, 
the radius of the dense spike can be large enough to accelerate the spaceship close to the speed of 
light. 

PACS numbers: 95.35.+d, 45.40.Gj, 89.30.-g 
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It is difficult for human to reach the stars using current rocket technology. The energy sources range from chemical 
fuel, nuclear power and even anti-matter conceptually The major problem in these systems is that propulsion required 
large amount of time and fuel[l[. We try to solve this problem by starting at the requirement of large amount of 
_ fuel. We all know that current rockets in function are chemical rockets which take oxidant and fuel at the same time. 
. Interestingly, the airplanes with similar propulsion system only take fuel, without oxidant, because in the atmosphere 
^\ ' there are enough oxygen which are absorbed by airplane engines during the flight. Similarly, if there are enough fuel 
Cn| . in the universe, the spaceship may absorb them during its flight like airplanes absorb the oxygen. Fortunately, dark 
matter is widely spread in the universe and the mass density is about five times of the baryonic matter density 
\ which make it a possible new energy source for interstellar flight. Thus the requirement of fuel may be solved in the 
00 ' self-help way with dark matter as the energy source. 

o 

II. DM ENGINE AND ACCELERATION IN THE SATURATION DENSITY 

We give an example of DM engine which uses DM annihilation remnants as propulsion. FigQ]is a sketch of the DM 
engine for this kind of new spaceship. The DM engine is the box in the picture. Here we assume the DM particle and 
the annihilation products can not pass through the wall of the box. In picture A, the space ship moves very fast from 
right to left. The DM particles, which are assumed to be static, go into the box and are absorbed in the picture B. 
In the picture C, we compress the box and raise the number density of the DM for annihilation, where we assume the 
annihilation process happens immediately In the picture D, only the wall on the right side is open. The annihilation 
products, for example Standard Model (SM) particles, are all going to the right direction. The processes from A to 
D are the working cycle for the engine. Thus, the spaceship is boosted by the recoil of these SM particles. Note the 
spaceship can decelerate by the same system when it reaches the destination, by opening the left wall in the picture 
D. 

This kind of new spaceship has a very interesting character that the faster it is, the easier it accelerates. In the 
picture A, we assume the rest mass of the spaceship is M and its velocity is (3 in the unit of speed of light. The time 
for one cycle of the engine is dt and the area of the engine is S. During one work cycle, the number of DM particles 
collected by the engine is N = (3dt ■ S ■ , where prj and m d are the density of DM and mass of DM, respectively. In 
picture D, we assume there are only one kind of particles X as the annihilation products for simplicity. The annihilation 
process is DD — > XX, with the mass mx < rrtrj. For the DM particles, we assume DM mass mj ~ O(10QGeV) and 
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FIG. 1: The illustration of work cycle for the DM engine. 



the annihilation products to be SM fcrmions mainly, which arc quite natural in SuperSymmctry and Extra Dimension 
models. Thus, the mass of the annihilation products mx arc quite small comparing with the mass of dark matter 
mo- So it is reasonable to use the approximation mx = 0, where products are treated as massless photon in the 
following calculation. Using the conservation of energy and momentum, we can get 



/ 0dt ■ S ■ p D + P° = P° + dp + e m 
\ 0+p 1 = p 1 +dp 1 -e-9 [) 

where p° = M7 and p 1 = M~/f3 are the energy and the momentum of the spaceship, and 7 = (1 — /3 2 )" 1 / 2 , e is 
the energy of the massless photons. 9 is defined as the propulsion efficiency which means 9 G [0,1]. For example, if 
the annihilation particles all go to the right direction, then 9 = 1. However, if the annihilation particles have equal 
possibility go into any direction in the right hemisphere, then 9 = 1/2. Moreover, it can also be used to count in the 
other inefficiency of the engine. By the Eqn[Tl one can get the differential equation for the velocity, 



k/3 3 df3 



(2) 



9- 1 + 13 ' dt 

where k = Spu/M. In the non-relativistic region, the above equation has the simple form, 



9k.p=f t . (3) 

To carry out the numerical calculations, we give some reasonable parameters first. We assume the weight of the 
spaceship is M = lOOton and the area is S — 100m 2 , according to the current rockets and space shuttles. For the 
DM density, we use the saturation density in the center of cusped halos p sa t- The saturation density in the halo is 
due to the balance between the annihilation rate of the DM [(crv) Psat/nriD]' 1 and the gravitational in falling rate 
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of the DM (Gp)- 1 / 2 , where the p is taken to be 200 times of the critical density. Thus the saturation density is 
p sa t ~ 10 19 Mq ■ fcpc _3 [3, Q- The propulsion efficiency is taken to be 9 — 0.5, since in the picture D we assume the 
annihilation particles have equal possibility go into any direction in the right hemisphere. We can also rewrite the 
parameter k as following, 



k = 2 x 10 _4 s 



4„-l 



S 100fon\ 



1O 19 M • kpc~ 3 100m 2 M ) 
One can solve the Eqn[2] and get the time and length needed for acceleration as the function of velocity, 
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where /3o is the initial velocity at the t — 0. We give the plots of the above equations in Fig[2] We can see the 
velocity increases exponentially with time, since the acceleration is proportional to the velocity. In the non-relativistic 
region, where (3 <C 1, the Eqnj5]and the Eqn[(5]can be simplified as 



= /3oe 



8kt 



(7) 



L=(0k)- 1 (p-p o )- 



(8) 



The initial velocity /3o is taken to be 10 _6 c which is much smaller than the first cosmic velocity. However, the result 
is not sensitive to the initial velocity, because of the exponential increase. In Figj2l we see the spaceship can reach the 
relativistic speed in about 2 days and the length needed for acceleration is about 10 -4 pc. From the above equations, 
if the DM density p and the area of the spaceship S are larger, the time t and length L needed for acceleration will go 
down. If the mass of spaceship M is larger, the time t and length L needed for acceleration will increase. However, 
the mass of DM particle does not have great influence on results, but the DM density does. 
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FIG. 2: The velocity as a function of the time and length needed for the acceleration in the saturation density. 
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III. ACCELERATION IN THE HALO OR SUBHALO 



Before celebration for the reach of relativistic speed, we should check whether the saturation region in the halo or sub- 
halo is large enough for the above calculation. The DM profile can be parameterized as p(r) = ( r / r yfw^jfk )q](/3-t)/o i 
where p s and r s are the scale density and scale radius parameters respectively. The parameters (a,/3, 7) are (1,3, 1) 
for NFW profile @. Since we are interested in the central region of halo, where r <C r s , the profile can be simplified 

as, 



This profile is singular at the center of the halo, ft is natural to have cut-off for this singularity due to the balance 
between the annihilation rate of the DM and the gravitational in falling rate of the DM. The saturation DM density 
is taken to be p sa t ~ 10 19 Mq ■ kpc~ 3 , thus the radius of saturation is r sat = p s r s / p sat . Once we know the scale 
density p s and scale radius r s , we can calculate the saturation radius r sat . The p s and r s can be fully determined by 
the concentration model and DM halo mass, which will be calculated in the appendix. Here we show the saturation 
radius r sat in Fig(3j We can see the saturation radius of halo or subhalo is much smaller than the required length for 
acceleration to the relativistic speed. 




The mass of (sub)halo (Msun) 



FIG. 3: The saturation radius r sat for different (sub)halo mass and concentration models. The B01 and ENS01 stand for 
different concentration models which are described in the appendix. 

In FigHJ we show the details of acceleration in the subhalo. The subhalo with mass 10 6 M Q in B01 model is taken 
as an example, which has saturation radius of about 10~ 9 pc. Starting from the center of subhalo with initial velocity 
(3q = 10 _6 c, it reaches the velocity of about 10 _5 c when it leaves the saturation region, which can be read out from 
the Fig|2] However, the rest of the subhalo is not sufficient to accelerate the spaceship to the relativistic speed, since 
the density begins to decrease by r . By solving the differential equation numerically, we can get the relations among 
velocity, time and distance in FigHJ We can see the spaceship reaches the velocity 10~ 4 c in about ten days. However, 
its velocity hardly increases after that, since the DM density goes down quickly. We can see that the acceleration is 
fastest in the saturation area of the halo. But the rest region of subhalo can still accelerate the spaceship from the 
velocity 10 _5 c to the velocity 10 _4 c. 

To better understand the acceleration power of the (sub)halos, we give the velocity at different times for different 
(sub)halos mass and spaceship parameters in FigjS] From the picture on the left, we can find that the subhalos have 
the power to accelerate the spaceship to velocity 10 _5 c ~ 10 _3 c with reasonable parameters S/M = 100m 2 /100<on. In 
the first few hours, the spaceship flics in the saturation region where they will be accelerated to velocity 10 _6 c ~ 10~ 4 c, 
which can be understood with the help of Fig|2]and Fig[3J Out of the saturation region, the velocity of the spaceship 
can get further boosted by about one order in the r _1 density region. Note that the above accelerations take place 
at the very center of halo, which is far less than the scale radius r s . The above results rely on the parameters of 
spaceship, e.g. the ratio S/M. If we lower the the weight of spaceship and increase the area of the engine, the velocity 
we can achieve will significantly increase. We specially give the plot on the right for S/M which is ten times larger, 
although the parameters maybe unreasonable in practice. It shows the the corresponding velocity increases about ten 
times. The main reason is the velocity at r = r sat increases by ten times, which can be understood with the EqnJ5] 
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FIG. 4: The velocity as a function of the time and length needed for the acceleration in the DM (sub)halo. The spaceship 
starts from the (sub)halo center with initial velocity f3o = 10~ 6 c. 
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FIG. 5: The velocity at time t = lday and t — lmonth for different (sub)halos mass and spaceship parameters. The concen- 
tration model is taken to be B01. The spaceship is still assumed to be started at the center of (sub)halo with initial velocity 
I3 Q = 10 _6 c. 

Anyway, the (sub)halos seem difficult to boost the spaceship to relativistic velocity because their saturation radius 
is small comparing with the required acceleration length 10 -4 pc in the FigO In the above calculation, we assume 
there are no baryonic matter in the halo. The gravity from the DM halo have negligible effects on the spaceship, even 
at the saturation region. Note the saturation density p sa t ~ 10 19 M Q • kpc~ 3 is much smaller than the density of water 
Ig/cm? ~ 10 31 M© ■ kpc~ 3 . 

However, in case there are baryonic matter in the halo, it may modify the DM profile. The adiabatic contraction 
due to dissipating baryons can steepen the DM profile Q. A more interesting case is that there is a central black hole 
in the DM halo, for example at galactic center. The DM density can become a dense spike due to accretion by the 
black hole, assuming adiabatic growth of the black holeQ- The annihilations in the inner regions of the spike set a 
maximal dark matter density p cor e = ( <T "° l)h ~ 10 17 M Q • fcpc -3 , where mo is the mass of DM particle, and tbh is the 
age of black hole, conservatively 10 10 yr. And more importantly, the radius of the core can be as large as O(10~ 2 pc) 
for inner cusped model like NFW profile. Recall the EqnlHl the required acceleration length for velocity 0.9c is about 
O(10~ 2 pc) in this case, which means the spaceship can achieve the velocity close to the speed of light. 
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IV. CONCLUSION AND DISCUSSIONS 



In this work, we give an example of DM engine using DM annihilation products as propulsion. The acceleration is 
proportional to the velocity, which makes the velocity increase exponentially with time in the non-relativistic region. 
The important points for the acceleration are how dense is the DM density and how large is the saturation region. 
The parameters of the spaceship also have great influence on the results. For example, the velocity will increase 
if S/M increases. We show that the (sub)halos can accelerate the spaceship to velocity 10 _5 c ~ 10 _3 c under the 
reasonable parameters of spaceship. Moreover, in case there is a central black hole in the halo, like galactic center, 
the core radius of DM can be large enough to accelerate the spaceship close to the speed of light. 

We have used three assumptions in this work. First, we have assumed static DM for simplicity. But the DM 
particle may have velocity as large as O(10 _3 c). Once we know the velocity distribution of DM, it can be solved by 
programming the direction of the spaceship when speed is low. An analogue in our daily life is airplanes work well 
in both headwind and tailwind. Second, we have assumed the DM particles and the annihilation products can not 
pass through the wall of the engine. For the annihilation products, they may be SM fermions which have electric 
charges. Thus we can make them go into certain direction by the electromagnetic force. The most serious problem 
comes from DM which are weakly interacting with matter. Current direct searches of DM have given stringent bound 
on cross-section of DM and matter. It may be difficult using matter to build the containers for the DM, because 
the cross-section is very small. However, the dark sector may be as complex as our baryon world, for example the 
mirror world. Thus the material from dark sector may build the container, since the interactions between particles 
in dark sector can be large. Third, the annihilation process is assumed to happen immediately in the picture C. This 

2 

is the second serious problem we should pay attention to. The annihilation speed takes the form, A = (av) £™£ ■ 

The (av) is taken to be the natural scale of the correct thermal relic, which is 3 x 10 _26 cm 3 s _1 . One can show that 
A = 2.2 x 10~ 7 cm~ 3 s~ 1 . However, the number density of the dark matter is no = = 4 x 10 9 cm~ 3 . Thus, 
to make the annihilation process efficient, we have to compress the volume of the engine to raise the annihilation 
speed. Whether it can be achieved in the future is not clear. Nevertheless, the engine works in the vacuum where the 
baryonic matter is dilute, which means we do not need to worry about the pressure from the baryonic matter. 

Sometimes, when looking at the N-body simulation pictures of DM, I think it may describe the future human 
transportation in some sense. In the picture, there are bright big points which stand for large dense halos, and the 
dim small points for small sparse halos. Interestingly, these halos have some common features with the cities on 
the Earth. The dense halos can accelerate the spaceship to higher speed which make it the important nodes for the 
transportation. However, the sparse halos can not accelerate the spaceship to very high speed, so the spaceship there 
would better go to the nearby dense halo to get higher speed if its destination is quite far from the sparse halos. 
Similarly, if we want to take international flight, we should go to the nearby big cities. The small cities usually only 
have flights to the nearby big cities, but no international flights. Thus we can understand the dense halos may be 
very important nodes in the future transportation, like the big cities on the Earth. 



APPENDIX: DARK MATTER HALO AND SUBHALO PROFILES 



Based on N-body simulations, the DM distribution can usually be parameterized as, 

P(r) = (r/r s )-<[l + (r/r s ) a ]^)/^ (AJ) 

where p s and r s are the scale density and the scale radius parameters respectively. The parameters (a, /3, 7) are 
(1,3,1) for NFW profile. In this appendix, we briefly introduce how r s and p s are calculated. The two parameters 
can be determined once we know the (sub)halo mass M v and the concentration parameter c„ which depends on the 
specific concentration model. The calculations are following the method in Ref.0]. In the appendix of Ref.0], we 
have shown how to determine the r s in detail. 
For the NFW profile, the scale radius r s is 



Cv(M v ) 

where r v is the virial radius of the subhalo which is often approximated as the radius within which the average 
density is greater, by a specific factor A = 200, than the critical density of the Universe p c = 139Mokpc~ 3 (Mq 

1/3 

is mass of the Sun). Thus r v can be expressed as r v = ( j^j§y££ J ■ The c v is the concentration parameter of 
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the subhalo which is determined by the subhalo mass M v and concentration model. We use the same method as 
Ref. [ij which adopts two concentration models, which are ENSOlQ and B0l[l(|- In the Ref. [H, the c v is fitted in a 
polynomial form as 



ln(co) = /,Cj 



i=0 



In 



My 



(A.3) 



where d = {3.14, -0.018, -4.06 x 10" 4 , 0, 0} and {4.34, -0.0384, -3.91 x 10" 4 , -2.2 x 10~ 6 , -5.5 x 10" 7 } for ENS01 
and B01 model respectively. 

The density scale p s can be determined by the mass relation J p s (r)dV — M v . One can get the scale density p s , 



% fw = M v /[4irr 3 s A(c v )\, 



(A.4) 



where A(c v ) = In(l + c„) — c„/(l + c v ). In Fig[6j the scale radius r s and the scale density p s are plotted as a 
function of subhalo mass M v to show how large and how dense the subhalos are. We can see the scale radius r s is 
quite large which shows the acceleration is mostly done in the r <^ r s region of the (sub)halo. 
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FIG. 6: The scale radius r 3 and the scale density p s as a function of (sub)halo mass M v . This plot assumes the (sub)halos 
have NFW profile. 
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